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EXECUTIVE SUMMARY 


Most studies on the effectiveness of interventions for reducing children's blood lead 
levels (PbB) have not distinguished declines in PoB due to program effectiveness from 
seasonal and age-related fluctuations in PbB. In this report, seasonal fluctuations and age 
effects in 1990-94 blood lead levels for a northern urban environment are studied, using 
data from 13,476 children screened for blood lead in Milwaukee, Wisconsin. The purpose 
was to determine whether there were seasonal and age trends, and if so, to estimate the 
magnitude of the trends. These estimates can then be used to help interpret studies with 
blood lead monitoring data, especially studies on the effectiveness of interventions to 
reduce blood lead levels. 


The Milwaukee data showed sizeable seasonal and age trends in Milwaukee children's 
PbB levels. Blood lead levels were about 40% higher in the summer than the winter, and 
about 15-20% higher at ages two to three years than at ages less than one year or ages 
five to seven years. Statistical methodology was developed to account for these 
fluctuations, so that the effectiveness of intervention programs may be quantified. These 
estimates are being used in studies of the effectiveness of lead interventions in Milwaukee. 
The methodology was described in considerable technical detail to facilitate analyses of 
seasonal and age effects in PbB in other environments. 


A tentative result suggests the magnitude of seasonal PbB fluctuations may be greatest 
for children less than four years old. Better understanding of the reasons for the trends 
might help to better determine mechanisms for reducing childhood lead exposure. 
Seasonal fluctuations in PbB are probably greater in cooler environments such as 
Milwaukee's, where seasonal changes in exposure to outdoor lead sources and sunlight 
are more extreme. At least in the northern U.S., the magnitude of the seasonal and age 
trends are large enough so that they must be considered in the design and interpretation 
of any blood lead monitoring results. 


xi 


1 INTRODUCTION 


Many Federal, State, and local programs have been implemented to reduce lead 
exposure in children. However, most of the studies of the effectiveness of the 
interventions for reducing blood lead (PbB) levels (see Burgoon, et al, 1994 and U.S.EPA, 
1995a) have ignored effects of seasonal fluctuations in PoB levels and dependencies of 
PbB levels on age. A few studies on abatement effects eliminated the confounding of 
seasonal effects but not age effects, by spacing PobB measurements one year apart. For 
retrospective studies, the requirement of many repeat measurements one year apart is not 
feasible. Thus, quantification of the effects of abatements and other promising intervention 
strategies has often not been possible. For example, a study on the effect of educational 
and counselling interventions for reducing PbB levels in children (Kimbrough, et al, 1994) 
reported that "educating parents proved a very effective tool" for reducing PbB levels, but 
did not estimate the decline in PbB levels due to the intervention. Instead, declines in the 
Granite City, Illinois children's PbB levels following interventions were simply reported as 
being too large to be attributed entirely to seasonal and age effects. 


In this report seasonal fluctuations and age effects in 1990-93 PbB levels for a northern 
urban environment are studied using data from 13,476 children screened for blood lead 
in Milwaukee, Wisconsin. The purpose is to determine whether seasonal and age effects 
exist, and if so, to develop simple adjustments to allow for quantification of the effects of 
education and abatement programs. Some previous analyses of data from the 1970's and 
early 1980's, such as U.S.EPA(1995b), indicated seasonal trends; this analysis sought to 
verify and estimate the current trend. Statistical methodology is described in detail to 
facilitate analyses of seasonal and age effects in PbB for other environments. Tables with 
adjustment factors specific to the 1990-93 Milwaukee PbB data are also included. The 
adjustment factors are necessary for planned analyses of the effectiveness of abatement 
and lead educational programs in Milwaukee from 1990-93. 


1.1 PEER REVIEW COMMENTS 


This study was reviewed independently by members of a peer review panel. 
Comments which are important for interpreting the study results or which had an important 
impact on the report are discussed below. 


Some reviewers wanted a clearer description of the Milwaukee Health Department 
screening program expansion, and how the expansion may have affected results of the 
analysis and validity of the seasonal adjustment factors. In response, portions of the text 
were rewritten. Table 1 was added; it delineates the changes and shows how the available 
data was used. Also Figure 1 now includes a needle plot showing the number of children 
screened for blood lead from 1983-1993. Geographical information was requested, so the 
appendices now include Table C1 which shows blood lead level results by zip code. 


One of the reviewers was concerned about correlations among residuals from different 
time periods. In response, greater emphasis was placed on Figure 12, a residual plot for 
90th percentile results which indicated the residuals were not correlated. The analysis of 
90th percentiles better describe blood lead levels of children at high risk for lead poisoning 
than the analysis of mean blood lead levels which indicated correlation among residuals. 
More discussion about this correlation was added to the text. 


One of the reviewers thought that the use of the Beta model was not well justified. In 
response, figures were added to show that the fit to the model was very good for the 90th 
percentiles. The beta model, with a minimum of parameters fit the data as well as 
sinusoidal models. The Beta model also allowed for direct assessment of features such 
as the potential for abrupt or asymmetric seasonal changes in PbB between winter and 
summer. A test for symmetry is now detailed in Section 4.2.1. 


Some of the reviewers were concerned about limitations inherent in making 
adjustments for removal of seasonal trends. Comments were made that well-designed 
randomized trials would likely be a better approach for studies of intervention 
effectiveness. In response, text in the Discussion section now more completely discusses 
these limitations, and mentions that some seasonal effects might be controlled for in some 
well-designed studies. However, formulating well-designed studies for evaluating 
effectiveness of interventions is problematic. Furthermore, retrospective studies, which 
by definition can not control for seasonality, have certain advantages. First, retrospective 
studies usually would include a larger number of children. Second, retrospective studies 
would not create artificial circumstances which could lead to invalid conclusions. Finally, 
they would not arbitrarily deprive control group children of benefits from interventions 
being evaluated. 


One of the reviewers was concerned about inferences about the interaction between 
age and other effects since the participation in the screening program may be different for 
older children (ages greater than three years). Although these are valid concerns, older 
children were usually screened for the same reasons as younger children, and only rarely 
because of clinical indications. Also, Figure 4 indicates that the relationship between age 
and PbB seems consistent from ages 2 to 7. 


2 DATA DESCRIPTION 
2.1 DATA COLLECTION 


The PbB data is a result of widespread blood lead screening of Milwaukee children, 
generally less than seven years old. The screening, occurring in many locations, attempts 
to identify children with elevated blood lead levels, so steps can be taken to reduce lead 
exposure and lead-related health impacts. The program began in the 1980's with a few 
health providers and laboratories reporting blood lead measurements to the Milwaukee 
Health Department (MHD). In late 1989, the MHD improved its computerization of its 
records. The screening program expanded dramatically in late 1991. From 1992 through 
1993 baseline measurements on 10136 children (about 5000 per year) were sent to the 
MHD. This corresponds to a coverage of about 50%, since about 10,000 children are born 
in Milwaukee per year. 


At some of the blood screening locations, all of the results are reported to the local 
health department. Although by law all lead screening data is being reported to the local 
health department, some sites only reported elevated blood lead levels in the past. 
Samples for the sites reporting only elevated blood lead levels have been excluded, since 
this might bias the estimates of seasonality of results. Biases in estimates of average 
blood lead levels in the population of children in Milwaukee may also result from a 
procedural change in October, 1991, after which all blood analyses directly measured the 
lead levels. Prior to October, 1991, some children were screened using FEP (free 
erythrocyte protoporphyrin) blood analyses (instead of blood lead). Follow-up blood lead 
measurements were made for these children having an elevated FEP level. FEP 
measurements were not used for any of the analyses. 


Available data for the analyses includes PobB measurements on 25,665 children from 
1986 to March, 1994. All PbB measurements in the data set are the baseline 
measurements which were made before any intervention by the MHD. For children with 
multiple PobB measurements, the data for this seasonality analysis only includes the first 
measurement. Blood lead levels from children with prior FEP measurements were also 
excluded. The next three paragraphs detail the other exclusion criteria. 


By 1990, some sites reported all measurements to the MHD. However, other sites 
tended to report only elevated blood lead measurements to the MHD, and data on 9581 
children from these sites was excluded (see Appendix for further details). This exclusion 
was made to reduce the effect on estimates that would be due to changes in the reporting 
of measurements to the MHD. 


The analysis, tables, and figures (except Figure 1) excluded pre-1990 data on 2,051 
children and 1994 data on 557 children. The pre-1990 data was excluded because it 
represented blood lead levels for a small vaguely defined set of children who tended to 


have abnormally high blood lead levels. Data from 1994 was excluded because of 
concerns that it may not have been complete. 


Table 1 identifies the MHD measurements on 13,746 children that had been used for 
this report. Tables 2 and 3 are frequency distributions of these children by year of 
measurement and age. 


Table 1. | Use of Blood Lead Measurements from the Milwaukee Health Department 


Number of PbB measurements 
All recorded baseline 25,665 
measurements from 1986 to 
3/94 
Providers report all Figure 1 25,665 - 9,581 
measurements = 16084 


Measurements from All analyses, tables, and figures | 16084 - 2051 (pre90) - 557 
1990 through 1993 except figure 1, and exceptions | (early 94) 
below. = 13,476. 


Recorded ages from 6 months Table 3, Figure 4, Section 4.3 13,476-572 
to 7 years = 12,904. 


Table 2. First Time Participants in the Milwaukee Blood Lead Screening Program from 1990-1993 by 
Year of Measurement 


Number of Observations 


1990 1431 
1991 2.466 





1992 5.260 
1993 4319 
13.476 1000 





Table 3. _ First Time Participants in the Milwaukee Blood Screening Program from 1990-1993 by Age 
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2.2 SUMMARY STATISTICS 


Semi-monthly means (aggregated twice a month) and 90th percentiles of 
untransformed PbB measurements were used for graphing the data and as inputs for the 
formal analyses. The 90th percentiles were generally preferred for the purpose of 
quantifying effects of interventions for reducing lead exposure of children with higher blood 
lead levels. This is because the blood lead levels at the 90th percentile are more similar 
to those of children receiving lead interventions than those at the mean, median, or other 
measures of central tendancy. Since there are 24 semi-monthly periods per year, and four 
years of data, 96 values of means and 90th percentiles for time periods from 1990 through 
1993 were statistically analyzed for most analyses. For analyzing effects of age on 
seasonality results (described in sections 3.2.4 and 4.3), semi-monthly values were 
calculated for each of four age groups from 1990 through 1993 for a total of 4*96=384 
means and 90th percentiles. 


The measurements used for calculating the semi-monthly statistics may be thought of 
as samples from populations of baseline blood lead levels that would be reported to the 
MHD. The populations change as actual blood lead levels change, and as the screening 
program coverage becomes more complete. For asemi-monthly period, the actual sample 
of measurements depends on a number of factors that include which children are tested 
at clinics that report all measurements to the MHD, the timing of the measurements, 
measurement errors, and so on. 


From 1990 through 1993, the only obvious systematic change in the number of 
measurements reported to MHD occurred in late 1991 (see Figure 1). Nevertheless, 
subtle changes associated with the expansion of the screening program may have resulted 
in gradual systematic changes in the blood lead level summary statistics. These changes 
may have had a non-negligible effect on estimates of long term declines in blood lead 
levels. In contrast, the gradual changes would have had less effect on estimates 
characterizing the seasonality of blood lead levels, because seasonal levels differ 
substantially (by about 40%) within a relatively short period of time (between winter and 
summer). 


2.3 DESCRIPTIVE STATISTICS AND FIGURES 


As described in section 2.1, the analyses were based on PbB measurements made 
after 1989. Figure 1, shows a striking difference in the characteristics of the 1986 through 
1989 versus the post-1989 time series of mean PbB measurements. The excess variation 
in mean pre-1990 PbB levels is primarily due to the much smaller number of children 
screened before the 1990's. The large drop in observed PbB levels in late 1989 may 
partially be a consequence of the expansion of the screening program. It is likely that 
before 1990 when there was an even greater need to target children with the highest blood 
lead levels, the relatively few participating health care providers may have served children 
primarily in areas of the city where blood lead levels tended to be higher. By 1990, 
participating primary health providers were more numerous, and areas of the city with 
differing blood lead level characteristics may have been more evenly represented among 
the children being screened. Although useful for estimating seasonal and age-related 
trends, the data here, especially through 1991, is limited in its ability to determine long- 
term trends in blood lead levels. The most reliable source for determination of long-term 
trends is NHANES II and NHANES III (Pirkle, et al, 1994) which provided data from 1975- 
1978 and 1988-1991. These surveys found that blood lead levels had decreased 
substantially during the 1980's, but that a significant number of children still had blood lead 
contents at levels widely considered as unhealthy. 


Figure 2 is a plot of summary PbB levels from 1990-93 aggregated by semi-monthly 
period. Open circles and diamonds denote raw means and distribution-free 90th 
percentiles for each semi-monthly period. Solid boxes and triangles denote smoothed 


means and 90th percentiles. Smoothing reduces short-term fluctuations caused in part by 
variation due to the small number of children sampled each period. The smoothed values 
were weighted moving averages over time. For time period t, the smoothed value was 
equal to 30% of the raw value at time period t, plus 20% of the sum of the raw values at 
time period t-1 (preceding) and t+1 (following), plus 10% of the sum of raw values at time 
periods t-2 and t+2, plus 5% of the sum of the raw values at time periods t-3 and t+3. 


The smoothed curves in Figure 2 clearly show seasonal fluctuations in PbB 
measurements with a peak around July or August and minimum values occurring in the 
winter. In the summer, the peaks are easily identified, but for some winters, the PbB 
measurements seem almost constant. Especially for 1993, the plot indicates seasonal 
fluctuations may be asymmetric. That is, the rise to peak levels may be more gradual and 
over a longer period of time than the decline to the lowest winter levels. However, a long- 
term decline in lead levels would accentuate the steepness of the seasonal declines in the 
fall and winter. Apparent seasonal fluctuations could have also been confounded with 
changes associated with the expansion of the screening program. From late September 
to early October, 1991, the number of PoB measurements rose from 112 to 235. The 
simultaneous sharp drop in PbB levels may have been partially due to the inclusion of 
lower risk children into the screening population eligible for the study. 


Figure 3 illustrates results from functional form fits for both means and 90th percentiles. 
The functional forms were used to clarify issues such as the possible asymmetry of the 
data. Details about the functional form fit are provided in section 3.2.2 and Appendix A. 
From the fit, peak PbB levels most likely occur in August. 


Figure 4 presents plots of mean and 90th percentile PbB levels by age. PbB levels 
increased rapidly before the age of 2 years, and then declined gradually thereafter. 
Similar results were indicated in at least two other studies. In the Sydney Lead Study 
(Cooney, et al, 1989), blood lead levels increased from birth to 18 months and then 
declined for ages 18 to 48 months. In the Port Pirie Study (Baghurst, et al, 1992), PbB 
levels peaked at age two years. Figure 5, a plot of mean PbB by gender of child, shows 
no detectable difference in PoB between males and females. 


Figure 6 shows smoothed plots of 90th percentile PbB levels for two different age 
groups. The plots suggest substantial seasonality from 1990 to 1993 for children less than 
3 years old, but only in 1992 for the older children. The extent to which seasonality 
depends on age is uncertain, because of the substantial seasonality shown for all age 
groups in 1992. Plots in Figure 7 suggest similar seasonality in PbB for ages 6 months to 
1 year, 1-1.5 years, and 1.5-3 years in 1990, 1992, and 1993. In 1991, PbB 
measurements had an observed peak in late summer only for children less than one year 
old. In 1991, the early peak in PobB measurements for ages 1-3 years may be due in part 
to procedural changes that may have occurred in the blood lead screening program. 
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Figure 1 Mean Monthly Blood Lead Levels and Number of Children Screened Based on Milwaukee Health Department Data from Sites that 
Report All Measurements. 
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Figure 2} Semi-monthly Arithmetic Mean and 90th Percentile Blood Lead Levels with Smoothed Estimates 1990-1993. 
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Semi-monthly Arithmetic Mean and 90th Percentile Blood Lead Levels with Functional Form Fit 1990-1993. 


Figure 3 


10 

















6-9 9-12 12-15 15-18 18-21 21-24 24-27 27-30 30-33 33-36 36-42 42-48 48-54 54-60 60-56 66-72 


Age Category (in months) 


-—s—# Mean Blood Level +—a—a 90th Percentile Blood Level 





Figure 4 = Arithmetic Mean and 90th Percentile Blood Lead Levels for 1990-1993 by Age Category (n=12,904). 
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Figure 5 Semi-monthly Arithmetic Mean Blood Lead Levels for Males and Females. 
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Figure 6 Semi-monthly Smoothed 90th Percentile Blood Lead Levels for Two Age Groups. 
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Figure 7. Semi-monthly Smoothed 90th percentile Blood Lead Levels for Children Under Three Years of Age. 
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3 METHODOLOGY 


Empirical smoothing and model fitting approaches were used to characterize seasonal 
and long-term trends in the data. 


3.1 EMPIRICAL SMOOTHING APPROACH 


The empirical smoothing approach included 1) calculating arithmetic means and 
empirical percentiles ("raw means and percentiles") of PoB measurements for each half- 
month period, 2) smoothing the raw means and percentiles using weighted moving 
averages, and 3) plotting both the raw and smoothed statistics against time. 


Plots of the raw means and percentiles over time are often sufficient to depict seasonal 
and long-term trends of time series. Differences in raw means and percentiles of the PoB 
levels are approximately unbiased estimates of differences in baseline blood lead levels 
of children covered by the MHD screening program. Thus, the January, 1993 minus the 
January, 1992 90th percentile (sample) PoB measurements would on average be equal 
to the actual drop from January, 1992 to January, 1993 in the 90th percentile PbB levels 
of children covered by the program. Nevertheless, moving averages were calculated, 
because sampling variation and random short-term fluctuations can mask trends. In 
general, the moving average of a time series y,, is given by: 


(1) Si = ZppoW Yi i t= Pt, -.., N-p. 


Here, the weights w,, j=-p, -p+1, ..., P, add to 1, and pis the order of the moving average. 
The greater the order and the more similar the weights, the greater the degree of 
smoothing, and the greater the reduction in sampling variation and short-term fluctuations. 
The weights used here were wWo=.3, W,=.2, W=.1, W3=.05, and w,; = w,. Assuming 
independence, the sampling variance for the smoothed statistics would be only about 
19.5% (or 100% * .3°+2*(.27+.1°+.05°)) of the sampling variance of the raw statistics. The 
moving average of order 3 helped discern trends occurring over time intervals of a few 
months or longer. Graphs of the raw means and percentiles showed that PbB is higher in 
the summer than the winter, but may reach a peak in some years as early as April. Graphs 
of the smoothed means clearly showed that peak PbB's probably occur sometime in the 
summer, perhaps in July or August. 


In theory, the moving average may have also obscured observation of interesting 
sudden rises in blood lead levels occurring within a couple of months. In Figure 2, the 
moving averages rounded the observed rise in PbB from June to August, 1993, so that the 
August measurement appears as part of a more gradual rise in PbB from January to 
August. Although a sudden two-month 20% rise in PbB seems unlikely, the moving 
average would have obscured observation of the rise, if it did occur. In general, weights 
are chosen to balance the opposing objectives of 1) minimizing deviations between the 
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smoothed and raw data so that important shorter-term fluctuations are not obscured, and 
2) reducing unwanted short-term fluctuations to facilitate better discovery of long-term 
trends. Appendix A demonstrates, through a more mathematical description of the 
problem, how the chosen weights achieved this balance rather well. 


3.2 MODELLING 


The graphs of the smoothed statistics led to many questions about the size and timing 
of the seasonal fluctuations. Modelling was needed for testing hypotheses (such as 
whether the seasonality is symmetric), and to construct meaningful parameter estimates. 
Modelling was also used, almost as an exploratory tool as described in section 3.2.3, to 
investigate the complicated relationship between age and seasonality in PbB. The 
analyses used a novel nonlinear regression approach based on the beta function. Asa 
check, the data was also fit using a more common sinusoidal function model. Both the 
beta and sinusoidal functions models assumed that non-age related patterns in PbB could 
be reasonably expressed as the sum of an overall (downward) linear trend and a function 
for seasonality. 


3.2.1 Sinusoidal Functions 


The sinusoidal function model is given in equation 2, where summary PbB levels for 
the ith time period, Y,, are expressed as a linear function of sine and cosine functions: 


(2) Y, = & + Bot, + 2-11, o(B,jSin(2m1jx;) + B,,cos(2mjx;)) + e; 


Here, tis time in years ranging from t=0 on January 1, 1990 to t=4 on January 1, 1994, and 
x = the fractional portion of t so that x ranges from 0 on January 1 to 1 on December 31; 
e; represents the random error term. The e, would be independent if the e; represent 
differences between sampled summary PbB levels and "true" means or percentiles for 
populations of PbB levels of Milwaukee children. However, the random errors could be 
correlated if the e; also reflected the unpredictability of changing climatic conditions. For 
example, June and July PbB levels could be correlated if a hotter-than-normal July often 
follows a hotter-than-normal June, and PbB levels rise with temperature. 


The random errors might also be correlated if the model (in equation 2) is misspecified. 
For example, for the interval from 1990 through 1993, the long-term trend might be 
concave (turning downward). Then the random error terms as defined through equation 
2 would tend to be positive somewhere in the middle of the time interval, negative 
elsewhere, and the correlation would likely be positive. 


16 


SAS's PROC REG (SAS, 1990) was used to fit the data to the sinusoidal model, and 
higher frequency terms were eliminated using the appropriate F-statistic. The e, were 
assumed to be asymptotically normal, and the first-order correlation was calculated to test 
the assumption of independence of the successive e,. 


3.2.2 Beta Function 


The beta function was preferred for modelling fluctuations in PbB levels, because with 
a minimum of parameters, it allowed for features suggested by our empirical analysis. The 
seasonal component may not be symmetric, and it was noted that during winter months, 
lead levels seemed relatively constant. The beta function allows for this relatively constant 
low period with a minimum of parameters whereas the sinusoidal function does not. Unlike 
the sinusoidal form, the beta function assumes that the peak and minimum PbB level are 
reached once each year, and that between the times of peak and minimum PbB levels, 
(mean or 90th percentile) PoB levels change monotonically. The beta functional form for 
semi-monthly summary lead levels Y;, shown in equation 3 with subscripts omitted, 
includes components for both linear trend (L) and seasonality (S). 


(3) Y=L+S()+e where, 
L=a+ft, and 
S() = A(z/R)"® ((1-z)/(1-R))TOP, 


Here a, B, A, R, 6, and T are parameters, t=time in years, and z=the fractional portion of 
(t-), i.e. ||(t-d) - int(t-)||, so that if d>=0, z ranges from 0 on January 1 to 1 on December 
31. eis the random error component. @ is the phase parameter ranging from -0.5 to 0.5 
(-0.5 < @ < 0.5). The phase parameter is included to allow the model to fit the seasonal 
maximum at the appropriate time of year. The seasonal component equals 0 when z = 0, 
and reaches its maximum, A, when z = R. Thus A, (A > 0), is the difference between 
maximum and minimum values of the seasonal component. Note that when A=0, the 
model suggests there is no seasonal variation. Maximum lead levels occur on the 
365(R+)th day of the year. R, (0 < R < 1), determines whether the seasonal component 
is symmetric. If R=.5, the time between maximum and minimum lead levels is one-half 
year. If R > .5, the rise to maximum lead levels is more gradual than the decline to 
minimum levels. T (T > 0) determines the abruptness of changes in the seasonal 
component around the seasonal peak, the 365(R+) day of the year. Assuming that 
maximum lead levels occur during the summer, a large value of T indicates almost 
constant lead levels in the winter and most of the spring, a rapid rise to their peak in the 
summer, and then a rapid decline to the winter levels. A small value of T indicates less 
abrupt changes between peak and minimum lead levels. 


Parameters a, B, A, R, and T were estimated through nonlinear regression. The model 
was first fit assuming @=0 and the e, were independent. The model was then fit using a 
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range of values for o to test whether o=0 and to evaluate the sensitivity of inferences to 
changes in @. The independence assumption was checked by calculating a version of a 
correlogram from the studentized residuals. The correlogram is a plot of sample 
correlations g,, ... , Qo, where 


(4) g, =(2ryfy,)/(n-k) = t= 1, 2, ... 96-K+1, 


and the r, are the studentized residuals. Under the assumption of independence, the 
absolute value of the sample correlations would generally be less than 2/n° (Diggle, 1990). 


3.2.3 Accounting for the Expansion of the Program 

Two more nonlinear regressions were also considered to account for the greater 
coverage of the screening program in 1992 and 1993. First, weights were set equal to the 
semi-monthly sample sizes, to account for the possibility that the Var(e,) are approximately 
proportional to the sample sizes. The data was also reanalyzed after adding a term to the 


model to account for the effect of the program expansion in October, 1991. The term, 
denoted by P, is P = O if t < 1.79 (corresponding to October 1, 1991); = A otherwise. 


The sinusoidal model with the term for the program expansion was thus (omitting subscript 


ie 

(5) Y= + Bot + 2.1,5)(B,jSin(2m1x) + B,,cos(2njx)) + P +e 
The beta model with term P is: 

(6) Y=a+6t+P+S+e 


Weighted analyses (with weights equal to the sample sizes) based on the generalized 
models of equation 5 and 6 will be referred to as "weighted + P" analyses. 


3.2.4 Accounting for the Effect of Age 


As shown in Figure 6, both overall PbB levels and seasonality may depend on age. 
Equation 7 incorporates age effects into the model, so that 


(7) Y=y,(L+S,+P)+e where 
L=a+6t, and S, =A,(x/R)'™ ((1-x)/(1-R))'?) 


and y, and A, are parameters defined for age categories, 

k=0: (0.5,1], k=1: (1,1.5], k=2: (1.5,3], and k=3: (3,7) years. The age categories were 
chosen so that the age categories are similar with respect to: 1) the number of children 
within each category, and 2) the range of mean PbB levels associated with the ages within 
age category. From Figure 4, the range of mean PbB levels are about 2.5 ug/dL = (11.5-9) 
g/dl for the ages (0.5,1] years, 2 ug/dL = (13.5-11.5) g/dL for ages (1,1.5], 1.25 ug/dL 
= (14.75-13.5) g/dL for ages (1.5,3] years, and about 1.5 g/dL = (14.5-13) g/dL for ages 
(3,7] years. 


The y, are multiplicative age-related factors for overall PbB levels. For example, if y)=1 
and y,=1.3, children between 1 and 1.5 years would have average PbB levels 30% higher 
on January 1 than children up to 1 years old. After adjusting for overall lead levels, the A,, 
which represent seasonal differences between high and low PbB levels, are also allowed 
to depend on the same age categories. Results from this model (equation 7) are 
presented in Section 4.3. 


4 RESULTS 
4.1. BETA FUNCTION MODELLING RESULTS 


4.1.1 Unweighted Nonlinear Regression Analysis 


The minimum PbB occurred sometime between late October and early March, because 
the basic model of Equation 3 fit the data well (for details see Appendix A, Section A2) for 
o between -1/6 and 1/6. The date of maximum PbB, sometime between late July and mid- 
September, could be more precisely estimated, because lead levels apparently changed 
abruptly during the summer and early fall. For aminimum PbB occurring on January 1, the 
estimated maximum detrended PbB date (or the date of maximum PbB if the long-term time 
trend were removed) was August 13 with a 95% confidence interval (July 25, September 
2). For values of @ between -1/6 and 1/6, estimated maximum PbB dates ranged from 
August 4 to September 4. A, the difference in maximum and minimum lead levels, was 
also insensitive to choice of d. For o=0, the estimate of A from monthly mean PbB values 
was 3.66 with a 95% confidence interval (2.65, 4.67). For 90th percentile PbB values, the 
estimate was 6.70 with a confidence interval (5.22, 8.19). In 1993, the seasonal 
component would account for a 38% rise in mean PbB and a 40% rise in 90th percentile 
PbB levels from January to August. However, for some @, the seasonal component was 
symmetric whereas for other o it was not. Assuming a minimum PbB on January 1, the 
95% confidence interval for R would be .56 to .67, so that the seasonal component would 
be asymmetric. Results for other values of @ are shown in Table 4. Estimates of R ranged 
from .425 for a minimum PbB on March 1 to .88 for the minimum occurring on November 
1. 


Residuals for the unweighted analysis of means and 90th percentiles for @=0, are 
shown in Figures 8 and 9. Both figures suggest the residuals are generally larger for time 
periods before October, 1991 when sample sizes were smaller. Figure 8 also suggests 
a positive correlation in residuals for the semi-monthly means. Correlograms are shown 
in Figures 10 and 11. Correlations with absolute values above 2/(96)? = 0.204 are 
generally considered significant (see section 3.2.2). From Figure 10, the absolute 
correlations for semi-monthly means less than four months apart are significantly greater 
than zero. In contrast, it is not clear from Figure 11 whether the 90th percentiles are 
correlated. In Figure 11, one sample correlation of 0.23 at 3.5 months was (barely) 
significant (exceeded 0.204) and other sample correlations for short time lags (0.5 and 2 
months) were nearly significant. 


20 


Table 4. | Summary of PROC NLIN Results by Phase (0) 
MEANS DATA 


ee eee ee | 

| see | conte | tsa.0 | aaa | 

po Tans | se.07) | 671 | ccoptabte_| 
3/24 Acceptable 
4/24 (.37,.48) Acceptable 
sea | ante | |r 


90th PERCENTILE DATA 
Phase Date een aa Residual SS Fit? 


(82,94) 
(7.87) Acceptable 
(72,81) Acceptable 
(67.78) Acceptable 
Acceptable 
Acceptable 
(52,60) Accepiable 
(47.55) Acceptable 
(42,50) Acceptable 
(97.45) Acceptable 
(92,41) Acceptable 


MSE = 7.10 (for phase = 1/24) 





‘Maximum PbB occurs 365*R days after minimum PbB. R=0.5 implies seasonality is symmetric. 
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2 For means data, fit is acceptable (see Appendix A) if Residual SS < 276.1 = 264.9 + 3.81*2.92. 
3 For 90th percentile data, fit is acceptable if Residual SS < 666.1 = 639.0 + 3.81*7.10 
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Figure 8. Studentized Residuals from Unweighted Analyses of Means. 
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Figure 9. Studentized Residuals from Unweighted Analyses of 90th Percentiles. 
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Figure 10. Correlogram from Unweighted Analyses of Means. 
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Figure 11. Correlogram from Unweighted Analyses of 90th Percentiles. 


Correlation in residuals for means and 90th percentiles could be a result of correlations 
in short-term non-seasonal PbB fluctuations, and/or model misspecification. Short term 
non-seasonal fluctuations may include three components. The first component, sampling 
variation, would be approximately independent for different time periods. The second 
component is generated by unpredictable changes in climate and other factors that affect 
exposure and physiology and may cause short-term unpredictable changes in overall 
levels of PoB. Most of these unpredictable changes are short-lived, but some could last 
for several months, resulting in correlated means and 90th percentile PbB levels. The 
third component would include effects of short-term changes linked to the expansion of the 
PbB screening program; these changes may affect the types of children whose PbB 
measurements are reported to the MHD. Enough children would have to be sampled each 
month over a sufficiently extensive time period to be able to detect these correlations; 
otherwise the sampling variation would overwhelm the variation from the last two 
components. 
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Model misspecification might have resulted in correlated residuals. This may have 
occurred if from 1990-93 the long term trend was not strictly linear (see section 3.2.1). 
Nonlinearity in observed trends could be due to either nonlinear trends in summary values 
of PbB levels of all Milwaukee children, or nonlinear effects linked to the screening 
program's expansion. As discussed in the next section, the observed correlations seem 
partially attributable to the effects of an abrupt expansion of the MHD screening program 
in October, 1991. 
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4.1.2 Weights and the 1991 Procedural Change 


The results described in the previous section suggested that a weighted regression 
analysis would be appropriate, since the residuals tended to be larger before the program 
expanded in October, 1991. Weights were set equal to semi-monthly sample sizes, a 
proper choice if sampling error accounted for almost all of the random fluctuations. A 
weighted+P regression analysis was also performed to account for a possible October, 
1991 shift in the types of children included in the screening program. Residuals from the 
weighted+P fit to the 90th percentile data are shown in Figure 12. The plot shows no 
obvious pattern, indicating that under the model with the expansion term (See equation 6), 
the random error components for the 90th percentiles may be independent. Thus, the 
weighted+P analysis apparently yields valid confidence intervals for parameters R and A 
characterizing the seasonal variation of 90th percentile blood lead levels. 


Residuals for the weighted+P fit to the means data, shown in Figure 13, show a slight 
quadratic trend. The trend in the means residuals may be due to a nonlinear trend in 
mean PbB values of all Milwaukee children. Alternatively, nonlinear changes in the 
observed mean PbB values may be partially attributable to changes in the types of children 
covered by the program. It is not clear why the long-term trend seems linear for 90th 
percentiles, but may be nonlinear for the means data. It is possible that the difference in 
trends may be partially attributable to MHD interventions which target children at high risk 
of lead poisoning. 


Parameter estimates resulting from the unweighted regression, the weighted (without 
the expansion term), and the weighted + P regression analyses are shown in Table 5 
under the assumption that d=0. The first two columns of Table 5 compare parameter 
estimates from weighted and unweighted regression analyses. The choice of weights 
causes negligible to about 10% changes in estimates of slope, A (the difference between 
seasonal maximum and minimum PbB), and R (defines time of maximum PbB). The last 
column shows confidence intervals from weighted + P analyses. The estimate of A was 
substantially and significantly different from O only for the 90th percentile data. The 
similarity between the three sets of confidence intervals for A and R shows that the 
uncertainty about proper weighting of the data and treatment of the effects of 1991 
procedural changes may have had only a minimal effect on the results. 


Table 5. Comparison of Nonlinear Regression Results for @=0 


| MEANS DATA | 
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SS Le 
dali 
Slope’ -125,-.086 
poster) | coor) | c0t.674) | 
ts 65) 


90th PERCENTILE DATA 


Confidence Interval 
Parameter ; 
Unweighted Weighted Weighted + P' 


24.9278) (24.3270) (24.8272) 


Slope (-.153,-.114) (-.146,-.109) (-.105,-.050) 
Amplitude (5.21,8.19) (5.95,8.31) (5.25,7.49) 
(.564,.646) (.596,.654) (.602,.661) 

(2.65,9.64) (4.18,10.9) (4.06,10.4) 


ey Pe 


‘Weighted analysis using model with procedural term. 

2 

ug/dL 

3ug/dL per semi-monthly time period 

“Point estimate of slope for 90th % from weighted + P = -.0777 
Point estimate of procedural discontinuity in 90th % = -3.727 
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Figure 12. Studentized Residuals of 90th Percentiles from Weighted Analyses Using the Program 
Expansion Term. 
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Figure 13. Studentized Residuals of Means from Weighted Analysis Using the Program Expansion 
Term 
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1991 1992 
4+—1—« Sinusoidal 90th PCT e—s—s Beta 90th PCT 
-—s—# Sinusoldal Mean o>  BetaMean 





Figure 14. Predicted Values Derived from Fitting 90th Percentiles and Means to the Sinusoidal and 
Beta Models 
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4.2 SINUSOIDAL FUNCTION MODELLING RESULTS 

The model of equation 8 provided an adequate fit to the 90th percentile data. As was 
the case with the beta function, residuals resulting from the sinusoidal fit to the means data 
were somewhat correlated. 


(8) y,= a + Bot; + B,sin(2m1x;) + B,cos(2mx;) + P; + e, 


Higher frequency sine and cosine terms (see equation 2) contributed little to the overall 
fit. Parameter estimates from a weighted analysis are shown in Table 6. 


Table 6. | Parameter Estimates and First Order Correlation of Residuals for Weighted Sinusoidal Model 
Fit. 


goth Percentiles 
287 
-0786 

1 


ee ey 


‘ug/dL 
ug/dL per semi-monthly time period 





4.2.1 Comparison with Beta Function Modelling Results 


Table 7 compares sinusoidal and beta function modelling estimates of the 1) slope, 2) 
the difference between seasonal maximum and minimum PbB, 3) the date corresponding 
to 365*R, the date of maximum (detrended) PobB. The weighted (weights equal to the 
number of children) mean square errors (MSE), a statistic that indicates how closely the 
model fit the data, are also given. Plots of the predicted mean and 90th percentile blood 
lead values are shown in Figure 14. Both Table 7 and Figure 14 show similar results from 
the Beta and sinusoidal models. This may be an indication that the two models will yield 
reasonable results if 1) PbB levels change monotonically between peak and minimum 
levels, 2) seasonal fluctuations are not highly asymmetrical, and 3) changes in PbB around 
the peak level are not too abrupt. 


Although the choice of model (sinusoidal or beta) had only a minimal impact on these 
estimates, the beta functional model is a recommended tool for other PbB analyses 
because it can: 1) more directly assess asymmetry, and 2) more flexibly account for abrupt 
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seasonal changes. Table 8 illustrates how the beta model can be used to test whether the 
seasonal effect is symmetric, or if R=0.5. Parameter estimates and a weighted version 
(weights equal to the sample sizes) of the residual sum of squares (RSS) for R constrained 
to equal 0.5 (using the weighted+P approach for 90th percentiles) is given in the first 
column. The same are then shown for R unconstrained in the second column. Results of 
an F test, shown at the bottom of the table, indicate that there is insufficient evidence to 
reject the hypothesis that the seasonality is symmetric. Note that the parameter estimates 
in both columns of the table are virtually identical. 


Abruptness and asymmetry of seasonal fluctuations could depend on climate and other 
factors. It would be much more difficult to test the symmetry of the seasonality using the 
sinusoidal model. 


Table 7. | Comparison of Sinusoidal and Beta Function Modelling Results' 


MEANS DATA 


Sinusoidal Model Estimate Beta Function Model Estimate 


Date of Maximum PbB 
Weighted* MSE 


= 
° 
o 


pe) 
o 
oO 
pe) 
(op) 
oO 


90th PERCENTILE DATA 


0.078 


'Weighted+P analysis 

ug/dL per semi-monthly period 

3Absolute difference between maximum and minimum seasonal PbB levels (g/dL) 

“Weights equal to number of children; denominator equals 96-6=90 for sinusoidal model, 96-7=89 for Beta 
model 


fo) 
iN 
“I 
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Table 8. Testing Symmetry (R=0.5) through the Beta Model 


Results for 90th Percentiles 


Parameter 


A 
T 


N lo Jo [2 SS 
ne) ol (o>) oO 
Ni |o | Im nN 
B lo Jy [wo @ 
QD ID lo |N oO 


Weighted RSS 46830 


F-statistic = (46830-46760)/(46760/89) = 0.14 





4.3 AGE'S EFFECT ON SEASONALITY 


Results from this section will show that the seasonality and overall level of PoB may 
depend on age in a fairly complicated way. Parameter estimates for the model of equation 
6 (section 3.2.3) are shown in Table 9. Note from the intercept estimates that winter PbB 
levels increase with age, but the seasonal difference between highest and lowest PobB 
levels is greatest between ages 1 to 3 years. For children aged 1.5 to 3 years in 1993, 
there was an estimated 38% increase in PbB levels from January 1 to peak levels in 
August. This compares to a 30% increase for the youngest, a 41% increase for 1 to 1.5 
year olds, and only a 15% increase for children over 3. Estimated percent increases in the 
90th percentile PbB's were (from youngest to oldest) 30%, 48%, 37%, and 3%. The 
difference in the size of the overall seasonal trend from 1990 to 1993 by age group is 
striking. In fact, the observed seasonal trend for children over 3 years old was not 
significant. Summer PbB levels were actually higher for the 1 to 3 year olds than for older 
children. 


However, as Figures 15 and 16 show, conclusions about the dependency of the 
seasonality on age require caution. Only PobB measurements for ages 1.25 to 2.5 were 
consistently higher in the summer than the winter from 1990 through 1993. Similar 
"seasonal" patterns can be observed in 1991 to 1993 PbB measurements for ages 0.75 
to 1.25 years, and 1992-1993 PbB measurements for ages 0.5 to 0.75 years. Forages 2.5 
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to 4 years, "seasonal" patterns can be observed in 1992 to 1993, but not in 1990 or 1991. 
For ages 4 to 7 years, PbB was higher in the summer than the winter in 1992, perhaps in 
1993, but not in 1990 or 1991. 


The lack of any discernable pattern in PoB measurements for some age groups in 1990 
and 1991 may be largely due to sampling variation. (For children less than 9 months data 
was so limited the graph was omitted.) Conversely, observed patterns in the empirical 
data that suggest seasonality may be the result of short-term random fluctuations in PbB, 
sampling variation, or aberrations due to pre-1992 procedural changes affecting the 
reporting of PoB measurements. The inconsistent patterns in the PbB data for older 
children only demonstrate the need for analyses of data beyond 1993. The additional data 
would be helpful to determine whether patterns observed in older children in 1992 are 
evidence of real seasonality in PbB levels, or merely an aberration caused by random 
fluctuations and sampling variation. If the additional data would show substantially less 
seasonality for older children, adjustments described later in Section 4.4.1 could then be 
modified to account for the dependency of seasonality on age. 
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1990 1991 


6—o—s 0,75 to 1.25 Years (n=1842) 2.5 to 4.0 years (n=518) 
1.25 to 2.5 Years (n=821) 4.0 to 7.0 Years (n=431) 





Figure 15. Semi-monthly Smoothed 90th Percentile Blood Lead Levels by Age Group 
(1990 to 1991) 
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1992 1993 


-—e— 0,50 to 0.75 Years (n=399) ©—e—se 2,5 to 4,0 years (n=1444) 
-—s—+ 0,75 to 1,25 Years (n=3825) 6 4.0 to 7.0 Years (n=1873) 
1.25 to 2.5 Years (n=1631) 





Figure 16. Semi-monthly Smoothed 90th Percentile Blood Lead Levels by Age Group 
(1992 to 1993) 
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Table 9. Estimates from Model Incorporating Age 


MEANS DATA 


Age Group (in Years 


Parameter 
5, 1 1,1.5 1.5,3 


Intercept 15.2 16.5' 19.1 
PbB on 1/1/90 


Slope’ 


Seasonal 3.0 4.46 4.8° 
Difference* 


5.0 


1.09 


90th PERCENTILE DATA 
Age Group (in Years) 


Parameter 
(.5, 1] (1, 1.5] (1.5, 3] 
—— 
Slope 


Seasonal 5.9 
Difference 
.605 


-5.0 


Y 1.08 


‘Equals ay, = 15.2*1.09 

ug/dL per semi-monthly period 

SEquals By, = -.073*1.09 

“ug/dL 

SEquals A,*y, = 4.04*1.09 

°38% increase from January, 1993 levels 
"Fixed 
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4.4 ADJUSTMENT FACTORS 


Several methods for seasonal adjustment of data were considered for analyses of 
abatement and educational outreach programs. Adjustments were based on results for 
90th percentiles (instead of means), because the 90th percentiles better describe the PbB 
levels of children targeted by the intervention programs. Adjustments could be additive 
or multiplicative. For example, by defining additive seasonal adjustments s, for each of 24 
bimonthly time periods, the adjusted PbB, Y , of ameasurement Y taken during time period 
i would be Y = Y + s,. For multiplicative adjustments, Y = Y*s,. For our data, "average" 
PbB levels were only about twice as large during the summer of 1990 as compared to the 
winter of 1993, so either type of seasonal adjustment would yield similar results. Our 
model fitting procedures assumed that the effects of seasonality were additive. 


Multiplicative adjustments would have been indicated if the studentized residuals had 
been larger for higher predicted PobB measurements. The graph of the studentized 
residuals from the weighted + P analysis of 90th percentiles, shown in Figure 11, indicates 
no change in the size of residuals from 1990 to 1993, despite the decline in PoB. Thus, 
additive adjustments seem adequate. 


Adjustments could be calculated as simple differences in raw statistics, or based on 
model fitting results. Advantages of adjustments based on raw differences are their 
simplicity and "unbiasedness", but "raw" adjustments are often unstable. For 90th 
percentiles denoted by z, for the reference period and z, for period i, the simple adjustment 
would be §; = Z)-Z;. 


Adjustments based on model fitting results would be s; = Y, - Y,, where perhaps 
Y, = a + Bt, + Pp, + (6). 


Model-based adjustments would tend to be more stable, because the models are 
calibrated using all available data, as opposed to data from only the ith and reference 
periods. Model-based adjustments are recommended when data is limited, and the 
variance of raw estimates is very large. However, model fitting introduces an additional 
potential source of bias caused by model misspecification. Although for our data, results 
were almost identical for the beta and the sinusoidal models, adjustments would 
sometimes be highly dependent on the choice of model. 


Moving averages often offer a reasonable compromise between adjustment strategies 
based on model fitting and raw statistics. Seasonal adjustments shown in Appendix D are 
based on the moving average of order 3 described earlier in Section 3.1. Data used for 
the adjustment are within 3 periods of the adjusted period, yet the sampling variance may 
be less than 20% of the sampling variance for the raw summary statistics. Adjustments 
based on higher order moving averages would have been recommended if there had been 
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more sampling variation. Moving averages require the assumption that second derivatives 
(which describe the rates at which the slope of a function changes) of "true" seasonal and 
time trends are sufficiently small. 


4.4.1 Age and Seasonal Adjustments for Use in Studies of Intervention 
Effectiveness in Milwaukee 


Multiplicative age and additive seasonal adjustments used in a study of intervention 
effectiveness (U.S. EPA, 1996) in Milwaukee are given in Appendix D. For that study, the 
adjusted PbB values were "equivalent" PbB values for ages 1.75 to 2 years for 
measurements made in January, 1993. 


The adjusted PobB measurement, y , for a child in age group k during time period i with 
PbB measurement y was calculated using the formula: 


(9) y'=(y+A\)*M,. 


Here, M, is the multiplicative age adjustment factors (for the kth age group); A; is the 
additive seasonality adjustment factor (for the ith time period). The age and seasonal 
adjustment factors are shown in the last columns of Tables D.1 and Table D.2 respectively. 
For example, to calculate the adjusted measurement corresponding to a March 7, 1992 
measurement of 20ug/dL for a 1.4 year old child, note that the age adjustment factor (from 
Table D.1) is 1.083, and the seasonality adjustment factor (from Table D.2) is 1.42. The 
adjusted PbB would be (20+1.42)*1.083 = 23.2. 


Note that the seasonality adjustment precedes the age adjustment, and that the final 
adjustments depend to a small extent on the ordering. For the age adjustment to precede 
the seasonality adjustment, the seasonality adjustments would have had to be modified 
to represent the seasonal fluctuations for ages 1.75 to 2 years. The absolute difference 
in summer and winter PbB levels between ages 1.75 and 2 years might be greater than 
for other age groups. This is because seasonal differences seem proportional to average 
PbB levels, and PbB levels tend to peak around age 2 years. 


The multiplicative age adjustment factors are inversely proportional to the arithmetic 
mean PbB's based on data from 1990-1993. The additive seasonality adjustment factors 
are based upon moving averages of the "detrended" 90th percentile PbB's. Here, 
detrended means that the linear long-term trend had been removed from the semi-monthly 
90th percentiles before the moving averages had been calculated. The reason for this was 
to assure that the adjustments would only reflect changes in PbB due to seasonality. The 
"detrending" was designed to filter out other effects (such as effects related to the 
screening program expansion). The long-term trend was removed by first fitting the 
1990-93 90th percentiles to the Beta function model, equation 3, with d=0. From the 
model fit, the estimate of the downward trend was 0.1335yg/dL per semi-monthly time 
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period. Thus, the time series was detrended by adding 0.1335*i to the 90th percentile PbB 
measurements for time periods i = 1 ... 99. The time series was then smoothed using the 
moving average with weights .3, .2, .1, and .05. The additive seasonality adjustment 
factors, A,, were then set equal to the smoothed values minus the smoothed value for i=73 
(for the first half of January, 1993). 


It is difficult to determine a "best" method for filtering out all of the long-term (non- 
seasonal) effects on PbB from the seasonal adjustment factors. A control group should 
be used in studies of intervention effectiveness for reducing lead exposure, so that 
changes in adjusted blood lead levels due to the non-seasonal factors (that can not be 
filtered out) would affect both study and control groups. An alternative set of seasonal 
adjustment factors are givenin Table D.3. A desirable feature of these adjustment factors 
is the lack of any discontinuity due to the effect of procedural changes in October, 1991. 
For these adjustment factors, the time series of semi-monthly PbB values was detrended 
using the results from the weighted+P analysis given in the third column of table 3. The 
detrending was accomplished by first adding .07777%i to the 90th percentiles to remove the 
long-term trend. Then 3.727 was added for time periods after September, 1991 to remove 
the effect of procedural changes in October, 1991. Smoothing and the A, were then 
calculated as described in the previous paragraph. 
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5 DISCUSSION 


The Milwaukee data showed sizeable seasonal and age trends in Milwaukee children's 
PbB levels. Blood lead levels were about 40% higher in the summer than the winter, and 
about 15-20% higher at ages two years than at ages one year and five years. The 
seasonal fluctuations have been attributed to complex physiological changes linked to 
increased exposure to lead or increased sunlight in the summer. Exposure may increase 
in the summer because of factors that include increased outdoor playing time, more 
opening and closing of windows, increased hand-to-mouth activity, and drier leaded dust 
that more easily enters homes. In Milwaukee, the age effect, related to factors such as 
increased hand-to-mouth activity of two year olds, was similar to the age effect observed 
in at least two other studies, the Sydney Lead Study (Cooney, et al, 1989), and the Port 
Pirie Study (Baghurst, et al, 1992). 


Inconsistent results on seasonality in PbB from many small studies (see McCusker, 
1979) suggest that seasonal patterns in PbB differ by location and/or climate. The 
Milwaukee data shows substantial seasonal fluctuation in PbB, so these trends must be 
recognized in similar northern urban environments. Seasonal fluctuations in PbB are 
probably largest in cooler environments such as Milwaukee's where seasonal differences 
in outdoor play and exposure to sunlight are more extreme. The Milwaukee results 
contrast with suggestions of some researchers that early evidence of seasonal trends 
(Blanksma et al, 1969; Guinee, 1972) are no longer relevant because of the phaseout of 
leaded gasoline (see U.S.EPA, 1995a). 


The MHD data set has limitations that must be noted. The MHD data is routinely 
collected health department data, and was not subject to the types of data quality checks 
had the data been collected for other purposes. Also, the data is not reliable for estimating 
long-term trends before 1992, since blood lead screening was much more limited before 
1992. Thus, evaluating long-term trends in PoB was beyond the scope of this study (see 
NHANES III, Pirkle, et al, 1994). Nevertheless, the changes in PbB associated with the 
expansion of the screening program would have had less effect on estimates 
characterizing the seasonality of blood lead levels, because seasonal levels differ 
substantially within a short period of time (between summer and winter). 


Other studies should help refine our understanding of the factors that may cause 
complex seasonal and age-related patterns. Further study is suggested by a tentative 
result suggesting that the magnitude of seasonal PbB fluctuations may depend on age. 
For 1990, 1991, and 1993, the seasonal fluctuations in PbB levels, although substantial 
for ages 1-3 years, were limited for ages 4 to 7. 


A model based on the beta function, developed to analyze the complex seasonal 


patterns in PbB levels (see U.S.EPA, 1995b), should have broad application. Unlike 
traditional methods, the beta function model allows for direct assessment of features such 
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as the potential for abrupt or asymmetric seasonal changes in PbB between winter and 
summer levels. Statistical literature includes discussion of other time series with possible 
asymmetric rises and falls such as the Canadian Lynx data (Campbell and Walker, 1977) 
and sunspots data (Morris, 1977). Although evidence for asymmetry was limited in the 
Milwaukee data, it is necessary to be able to test for such features in PbB data from other 
environments, where seasonal characteristics in PoB levels may differ. 

Both empirical and model-based results indicated that seasonal and age related trends 
could influence results from studies of the effectiveness of interventions for lowering PbB 
in children. Although these trends might sometimes be controlled for in well-designed 
prospective studies, adjustment factors would be needed for trend removal in retrospective 
studies. However, retrospective studies have several advantages. First, they allow for 
adequate sample sizes, a major consideration because of the large variability in PbB 
levels. Second, they do not create artificial circumstances which could lead to invalid 
conclusions. They also do not arbitrarily deprive control group children of benefits from 
the interventions being evaluated. 


Seasonal and age-related adjustment factors were calculated through a four-stage 
process. First, 90th percentiles were calculated to summarize PbB levels for each of the 
96 semimonithly periods between 1990 and 1993. Second, the 90th percentile PoB values 
were fit to a model so that long-term trends in PbB could be removed. The seasonal 
adjustments were then based on the moving averages of the detrended 90th percentile 
PbB values. Finally, age adjustments were calculated as simple ratios of arithmetic means 
using predefined age categories. 90th percentiles were chosen as the summary statistics 
for the semi-monthly periods, because of their applicability to populations of children 
subject to PbB interventions. Since age and seasonal trends likely depend on many 
factors related to geographic location, type of environment (urban or rural), and time 
period, the adjustment factors shown in Appendix D are specific to Milwaukee from 1990 
to 1993. Nevertheless, the four-step procedure, with appropriate modifications, should be 
applicable to many other PbB data sets, so that the effects of abatement and educational 
interventions in other locations may also be quantified. 
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APPENDIX A. TECHNICAL DETAILS 
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A1. Inputs to PROC NLIN and Properties of Estimates 


SAS's PROC NLIN and the Newton-Raphson option were used for most' of the 
nonlinear regression beta function model fits. Details about the algorithm are provided in 
the SAS/STAT User's Guide, Volume 2. 


If the random errors e; are independent and identically distributed with finite variance 
o*, then the resulting parameter estimates of a, B, A, R, and T are consistent with 
asymptotic normal distributions. If in addition, the e; are also normally distributed, then 
upon satisfaction of certain regularity conditions, the estimates are maximum likelihood 
estimates. Details are provided in Seber and Wild (1989). 

Inputs into PROC NLIN include starting values for the parameter estimates, a model 
statement, bounds for some of the parameters, and first and second partial derivatives of 
E(Y) with respect to the parameters. The model, without a phase component, is: 

E(Y) = a + Bt + A(x/R)™ ((1-x)/(1-R))™®. 
The first derivative with respect to A is: 
d(E(Y))/dA = (x/R)'® ((1-x)/(1-R))™®. 


Let F = exp((In x - In R)TR + (In(1-x)-In(1-R))T(1-R)). Then the first derivative with respect 
to T is: 


d(E(Y)/dT) = AF(R In(x/R) + (1-R)In((1-x)/(1-R))). 


' A MATLAB program was written to generate the results shown in Table 8. 
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The other first and second derivatives can be gleaned from the programming code for 
the unweighted fit: 


proc nlin method=newton data=avlead; 
parms a0=17.3 a1=-.11 a2=5 a3=.6 a4=5; 
temp = (x/a3)**(a3*a4)*((1-x)/(1-a3))**(a4*(1-a3)); 
model y=a0 + (a1*t) + a2*temp; 

bounds a3>=0, a4>=0, a3<=1; 

der.a0 = 1; 

der.al =t; 

der.a2 = temp; 

temp2= log(x/(1-x))+log((1-a3)/a3); 

der.a3 = a2*temp*a4*temp2; 
temp3=a3"log(x/a3)+(1-a3)*log((1-x)/(1-a3)); 
der.a4 = a2*temp*temps; 

der.a0.a0 = 0; 

der.a0.a1 = 0; 

der.a0.a2 = 0; 

der.a0.a3 = 0; 

der.a0.a4 = 0; 

der.al.ai = 0; 

der.a1.a2 = 0; 

der.a1.a3 = 0; 

der.ai1.a4 = 0; 

der.a2.a2 = 0; 

der.a2.a3 = temp*a4*temp2; 

der.a2.a4 = temp*temp3; 
dfda3=temp*a4*temp2; 

der.a3.a3 = a2*a4*(temp/(a3*(a3-1))+dfda3*temp2); 
dfdt=temp*temp3; 

der.a3.a4 = a2*temp2*(temp+a4*dfat); 
der.a4.a4 = a2*temp3*dfdt; 

output out=preds p=yhat r=yresid; 
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A2. Incorporating Phase 
The model 
Y=a+ P(t) + S(o) +e 


was fit by comparing the residual sums of squares of repeated fits using a range of values 
for d. An approximate 95% confidence interval for @ would include all values of for 
which the residual sum of squares is no greater than RSS’ + 3.81*MSE.. Here, RSS’ and 
MSE’ are the residual sums of squares and mean square error values minimized with 
respect to $, and 3.81 is the critical value corresponding to a=.05 for the chi-square 
distribution with 1 df. 

Referring to table 5, RSS’ for 90th percentiles ~ 639.0 and MSE = 7.1 = 639/90, so the 
95% confidence interval for @ would include all values for which RSS < 666.1 = 639.0 + 
3.81*7.1. This approach allowed observation of how assumptions about © affected 
estimates of other parameters. 


A3. Properties of the Nonlinear Regression Estimates 


If the random errors e; are independent and identically distributed with finite variance 
o*, then the parameter estimates of a, B, A, R, T, P, are, upon satisfaction of regularity 
conditions (see Seber and Wild), consistent with asymptotic normal distributions. If in 
addition, the e,; are also normally distributed, then the nonlinear regression estimates are 
maximum likelihood estimates. 


A4. Mathematical Justification for Moving Average Weights 
The theoretical material in this section is from Diggle, 1990. Let [i(t) be a moving 


average for the time series y(t). Weights for a moving average might be chosen to 
minimize the quantity Q(a), where 


Q(a) = =X (y; - A(t))? + a f(H"(t))? at. 


The summation term, "the residual sum of squares", measures the closeness of the fit 
between the moving average and the original time series. The integral term measures the 
smoothness of the moving average. « determines the tradeoff between goals of obtaining 
a very smooth fit and minimizing the residual sum of squares. 


lf data are equally spaced at unit time intervals (as ours is), then Q(a) would be 
approximately minimized when 


(8) w, = .5*h"(K((i+.5)/h)+K((i-.5)/h)), where h=a?5, and 
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K(u) = .5 exp(-|u|/2") sin(.25r1 + |uj/2"). 
The weights of the 3-order moving average are within .05 of the weights given in (8) 


for a=1. For larger values of «, Q(a) would be minimized by a higher order moving 
average. 
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APPENDIX B. DATABASE DEVELOPMENT 


5A 


The Milwaukee Health Department (MHD) provided the data described in this report 
from the health department's lead case tracking system called "STELLAR". One of 
STELLAR's data files, LAB.BAS, contains records for each blood lead level reported to the 
health department. Each time a child's blood lead level was reported, all the children's 
blood lead levels in LAB.BAS were reviewed by MHD staff, and additional entries to the 
STELLAR data files were made when necessary. LAB.BAS also includes the 
corresponding date of measurement, the sample type, the child's name and date of birth, 
the child identifier, the address identifier, and the medical provider. As of July 1994, there 
were 75,084 records in this file. 


Medical providers and laboratories send data on children's blood lead levels to the 
MHD. Providers include primary care physicians, public health clinics, HEADSTART 
centers, and Women Infant and Children centers (WICs). Some providers reported all 
measurements and some only reported elevated levels. A list of all possible providers was 
created from the LAB file and each was called to verify the procedure used. All 
HEADSTART and WICs but only half of the clinics and physicians reported all 
measurements. Therefore, a final list of providers that reported all measurements was 
created and used in the creation of the database. The data analyzed for this report only 
includes measurements from providers who reported all measurements regardless of level. 


The database used for the analysis was created by first sorting the LAB file by the child 
identifier and sampling date. The first measurement in chronological order for each child 
was maintained in the database. Measurements were deleted if the corresponding 
provider name was not among those listed as reporting all measurements. The final 
database has 16,084 observations: 2,051 measurements occurring before January 1, 
1990, 13,476 measurements occurring between January 1, 1990 and December 31, 1993, 
and 557 measurements after December 31, 1993. 
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APPENDIX C. GEOGRAPHIC SUMMARY OF HEALTH DEPARTMENT DATA 
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Table C.1 summarizes blood lead data for each Milwaukee zip code which according 
to the 1990 U.S. Census had at least 2500 children less than 7 years old. U.S. Census 
population figures for each zip code are given in the first column. The table compares the 
Census data to the number of children who had blood lead measurements taken for first 
time (the number of children screeened). 


Estimates of coverage of the MHD screening program are given in the last column. 
The coverage of the MHD (screening) program for 1992-93 can be defined as the 
proportion of children who were are born in Milwaukee from 1992 through 1993 who have 
or will be tested for blood lead before their seventh birthday. The number of children born 
during 1992-93 would have been about 2/7 times the number of children less than seven 
years old in 1990. The number of children born in 1992-93 who would be screened by 
their seventh birthday would be approximately equal to the number of children tested for 
blood lead for the first time from 1992-93 under the following two assumptions. First, the 
number of children screened reached equilibrium by 1992. Second, the age distribution 
of the children screened was about the same for each year by 1992. A crude estimate of 
coverage for each zip code was calculated as the ratio of the the number of children 
screened during 1992-93 divided by 2/3 * (1990 Census population for ages up to 7 
years). 


57 


Table C1. Census Counts and Summary of Milwaukee Health Department Blood Lead Data by Zipcode’ 


Zipcode 1990 U.S Census Mean PbB Number Screened Coverage? 
Population g/dL 


Ages 0-7 Years 
1990-93 1990-91 1992-93 1992-93 


[saoos [zen | tao | ve | ot | som | 
[saos [roo | 50 | soa | 1100 | som | 
ae 
| scoop | ses | ta | ate | wes | tome 
| tot | erg | 


"Includes only zipcodes with >2500 children less than 7 years old. 
? Equals number screened during 1992-93 divided by (2/7)*population 
3 Over all zipcodes 
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APPENDIX D. 1990-93 MILWAUKEE PbB MEASUREMENT ADJUSTMENT 
FACTORS 
(These factors are not appropriate for other PbB data sets). 
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Table D1. Multiplicative Age Adjustment Factors for 1990-1993 Milwaukee PbB Data 


Age Category (Years Arithmetic Mean(ug/dL Adjustment Factor 





‘Includes 123 without any recorded birthdate, 106 with ages less than 6 months, and 237 with ages greater 
than 7 years. 
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Table D2. Seasonal Additive Adjustment Factors Without Procedural Correction for 1990-March, 1994 
Milwaukee PbB Data 


Year 1990 


Time 90th Detrended Smoothed Adjustment Adjustment 
Period Percentile 90th Detrended | Factors Without Factors* 
Starting PbB Percentiles’ Series? Detrending® 
Date 


Apr 16 0.482 
5276 
Aug 16 8.000 
2.684 


17958__| 17.943 7.045 0.503 
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Table D2 (continued) - Year 1991 


Time 90th Detrended Smoothed Adjustment Adjustment 
Period Percentile 90th Detrended Factors Without Factors* 
Starting PbB Percentiles’ Series? Detrending® 
Date 


2.344 
4276 
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Table D2 (continued) - Year 1992 


Time 90th Detrended Smoothed Adjustment Adjustment 
Period Percentile 90th Detrended Factors Without Factors* 
Starting PbB Percentiles’ Series? Detrending® 
Date 


1.252 
1.085 
4.797 
4,946 
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Table D2 (continued) - Year 1993 


Time Period 90th Detrended 90th Smoothed Adjustment Factors Adjustment 
Starting Percentile Percentiles’ Detrended Without Factors* 
Date PbB Series® Detrending® 


Jan 1, 93 18.000 17.440 0.000 0.000 


18.134 17.014 0.560 0.426 


1.509 


Apr 16 16.935 17.775 0.600 -0.335 


18.268 18.678 0.170 1.238 


May 16 21.202 19.642 -1.000 -2.202 


19.335 20.120 1.345 2.680 


Jun 16 20.469 20.709 -1.800 -3.269 


Jul 16 4.306 
Sop 16 4.980 
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Table D2 (continued) - Year 1994 


Time Period 90th Detrended Smoothed Adjustment Adjustment 
Starting Percentile 90th Detrended Factors Without Factors* 
Date PbB Percentiles’ Series? Detrending® 


‘Detrended by subtracting .1335*(73-i) from 90th percentiles. 
?Smoothed values of column 3 using weights .3, .2, .1, .05. 
3Adjustment factors if 90th percentiles had not been detrended. 
“Equals 17.44 (smoothed value at t=73) - column 3 values. 
5Uses 1994 90th percentile semi-monthly PbB values. 
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Table D3. Seasonal Adjustment Factors with Procedural Correction for 1990 - March, 1994 Milwaukee PbB 
Data 


Year 1990 


Time Period 90th Percentile Detrended 90th Smoothed Adjustment 
Starting Date PbB Percentiles’ Detrended Factors® 
Series® 
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Table D3 (continued) - Year 1991 


Time Period 90th Percentile Detrended 90th Smoothed Adjustment 
Starting Date PbB Percentiles’ Detrended Factors® 
Series® 


Jan 1, 91 21.370 20.555 0.612 


Jan 16 19.348 20.458 0.709 


21.726 20.401 0.766 


Feb 16 17.504 20.538 0.629 


22.781 21.331 0.164 


Mar 16 21.259 22.184 -1.017 


25.197 23,292 2.125 


Apr 16 21.814 24.004 -2.837 


24.092 25.147 3.980 


May 16 31.470 26.430 5.263 


24.047 26.097 4.990 


Jun 16 25.625 25.895 4.728 


27.203 25 893 4.728 


Jul 16 24.880 25.366 -4,199 


20.358 25,498 4.331 


Aug 16 28.336 26.252 -5.085 


27.814 26.418 5.251 


Sep 16 27.291 25.816 -4.649 


22.796 24.587 3.420 


Oct 16 23.774 23.470 -2.303 


23.551 22.340 4.173 


Nov 16 19.629 21.154 0.013 


19.707 20.522 0.645 


Deo 16 19.584 20.474 0.69 
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Table D3 (continued) - Year 1992 


Time Period 90th Percentile Detrended 90th Smoothed Adjustment 
Starting Date PbB Percentiles’ Detrended Factors® 
Series® 


Jan 1, 92 22.762 20.632 0.535 


Jan 16 18.940 20.530 0.637 


21.418 20.718 0.449 


Feb 16 20.095 20.820 0.347 


21.173 20,8683 0.904 


Mar 16 22.251 20.976 0.194 


20.328 20.898 0.929 


Apr 16 18.406 21.031 0.136 


May 16 2.634 


Jul 16 27.872 26.727 -5.560 


26.950 27.150 5,983 


Aug 16 28.028 27.178 -6.011 


27.108 26.985 5.788 


Sep 16 27.183 26.183 5.016 


24.261 25.165 3.994 


Oct 16 25.338 24.138 2.971 


21.416 22.966 4.798 


Nov 16 22.494 22.244 -1.077 


21.572 21772 0.605 


Dec 16 20.649 21.319 -0.152 
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Table D3 (continued) - Year 1993 


Time Period 90th Percentile Detrended 90th Smoothed Adjustment 
Starting Date PbB Percentiles’ Detrended Factors® 
Series® 


Jan 1, 93 21.727 21.167 0.000 


Jan 16 21.805 20.685 0.482 


19.282 19.982 1215 


Feb 16 18.960 19.490 1.677 


18.038 19.378 1.789 


Mar 16 20.116 19.846 1.322 


21.193 20.563 0.604 


Apr 16 20.271 aAdd 0.056 


21.549 21,959 0.792 


May 16 24.426 22.866 -1.699 


22.504 29.289 2.102 


Jun 16 23.582 23.822 2.655 


25.659 24.419 3.262 


Jul 16 24.037 24.747 -3.580 


23.816 25.275 4.108 


Aug 16 26.893 26.002 -4.836 


29.370 26.175 6.008 


Sep 16 25.248 25.198 4,031 


21.726 23.941 2774 


Oct 16 23.603 23.233 -2.066 


29,281 22.701 4.594 


Nov 16 20.859 22.069 -0.902 


21436 21.886 0719 


Dec 16 24.914 21.584 -0.417 
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Table D3 (continued) - Year 1994 


Time Period 90th Percentile Detrended 90th Smoothed Adjustment 
Starting Date PbB Percentiles’ Detrended Factors® 
Series® 





'‘Detrended by adding .0777*(i-73) and 3.727 for t>42. 
?Smoothed values of column 3 using weights .3, .2, .1, .05. 
SEquals 21.167 (smoothed value at i=73) - column 3 values. 
sizes (n=120). 
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APPENDIX E. 1990-96 MILWAUKEE PbB MEASUREMENT ADJUSTMENT 
FACTORS 
(These factors are not appropriate for other PbB data sets). 
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In response to one of the peer reviewer's comments a set of seasonality adjustments 
was calculated for log(PbB) values, and zipcode information for the screened children was 
incorporated into the trend analysis. During peer review, additional Milwaukee Health 
Department blood lead data became available through February, 1996. Seasonality trends 
were reevaluated with the new data, and adjustment factors were calculated for the 
extended time period January 1, 1990 through February 15, 1996. The basic methodology 
for calculating these new adjustment factors is as before. Methodological details are 
given in the following text. The new adjustment factors are given in Tables E1 though E3. 


Step-by-Step Details of the Adjustment Process 


Seasonal adjustment factors for the log transformed data were calculated through a six- 
step procedure. 


First, the log transformed data were adjusted to account for differences associated with 
the address zip codes of the screened children. This adjustment was done through an 
analysis of variance (ANOVA). The dependent variable was the log (PbB) value; the 
(three) independent classification variables were based on 1) the six bimonthly periods of 
the year (e.g., January - February, March - April, etc.), 2) (ten) half-year time intervals (first 
half of 1990, ..., last half of 1995), 3) zip codes (zip codes with less than 1000 children 
were combined). Each child's adjusted log(PbB) value was then set equal to the log(PbB) 
value minus the least squares mean value corresponding to the child's zip code. 


Second, the 90th percentile values were calculated for each semi-month period from 
January 1, 1990 through December 31, 1995. 


Third, the 90th percentile values were fit using the model: 


(11) log(Y;) = a + Byt; + P; + Si(o) + 
Bot? + By A, (t)) + BioAd(ti) +... + BigAg(ti). 


where t = time in years (for January 1, 1990, t=0). 


A,(t) = 4t for t < 0.25; A,(t) = 1 fort > 0.25. 
A,(t) = 2t for t < 0.5; A,(t) = 1 fort > 0.50. 
A,(t) = (4/3)t for t < 0.75; A,(t) = 1 fort > 0.75. 
A,(t) = t for t<1; A,(t) = 1 for t>1. 

A,(t) = 0.8t if t<1.25; A,(t) = 1 fort > 1.25. 
A,(t) = (2/3)t for t< 1.5; A,(t) = 1 fort > 1.5. 
A,(t) = (4/7)t for t<1.75; A,(t) = 1 fort > 1.75 
A,(t) = 0.5t if t<2; A,(t) = 1 fort > 2. 
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This is just an elaboration of the Beta model described in Section 3.2.2 and 3.2.3. Without 
the last five terms on the right-hand side, equation (11) is identical to equation (6) in 
Section 3.2.3. The quadratic term (8,t) was included to better account for nonlinearity in 
the long-term trend. The last four terms are included to account for medium-term 
nonseasonal fluctuations in the observed PbB values associated with the expansion of the 
screening program between January 1, 1990 and December 31, 1991. 


The fourth step is "detrending" the time series by subtracting out the components 
associated with the linear trend, P,, and the last five terms in equation 9. The fifth step is 
calculating a moving average of the detrended time series. Finally, the adjustment factors 
are set equal to the smoothed value at time t=73 (corresponding to January 1, 1993) minus 
the smoothed value for the time period when the child's blood sample was collected. 


The same procedure was used to calculate adjustment factors directly from the 
untransformed data, but with two exceptions. First, the PoB values were never adjusted 
using zip code information. Second (obviously), the log transformation was never used. 


The next paragraph describes the procedure for calculating age adjustment factors for 
the log transformed data. 


First, the log transformed data were adjusted to account for differences associated with 
the address zip codes of the screened children, and also differences associated with 
trends in PbB values over time. This adjustment was done through the analysis of 
variance (ANOVA) described above. Each child's adjusted log(PbB) value was first set 
equal to the log(PbB) value minus the least squares mean value corresponding to the 
child's zip code. Each zip code adjusted value was then adjusted for the time trend by 
subtracting the least squares value associated with the six month time interval (when the 
child's blood sample was collected). Average adjusted log(PbB) values were then 
calculated for thirteen predefined age categories. For each age category, the age 
adjustment was then set equal to the average for the age category minus the average for 
a reference category (1.75 through 2.0 years). 


Multiplicative age adjustments for unadjusted data were calculated by simply applying 
the exponential function to the adjustments for the log transformed data. 


As described in Section 4.4.1, untransformed data could be adjusted directly using the 
formula (equation 9, Section 4.4.1): 


y = (y+A)*M,. 
Here, M, is the multiplicative age adjustment factors shown in Table 3 (for the kth age 


group); A; is the additive seasonality adjustment factor from Table 2 (for the ith time 
period). 
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The seasonal and age adjusted log transformed data would be equal to the log 
transformed data plus the sum of the corresponding age (Table 3) and seasonality 
adjustments (Table 2). Examples for using the tables follow: 


Example 1: Suppose Pbb value on September 1, 1992 is 25. 
Then log(Pbb) = 3.219. 
log(Pbb) (adjusted for seasonality from Table 1) 
= 3.219-0.241 = 2.978. 
Seasonally adjusted Pbb = exp(2.978) = 19.6. (This represents the 
equivalent Pbb value for a measurement made January 1, 1993 using the 
adjustments that were based on the log transformed values). 


Example 2: Suppose Pbb value on September 1, 1992 is 25. 
Adjusted PbB = 25-5.729 = 19.3. 
(This represents the equivalent Pbb value for 
measurement made in Jan., 1993 based on the analysis of 
the untransformed values). 


Note that the adjustments from Tables E.1 and E.2 do not yield identical results. 


Adjustments from Table E.1 are being used for a retrospective analysis of paint 
abatements in Milwaukee. 
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Figure 17. Modelled (solid line) and Smoothed (dashed line) 90th Percentile Blood Lead Levels Using 
Log Transformed Data from 1990 to 1996 
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Figure 18. Seasonal Adjustment Factors (1990-96) for Log Transformed Data 


Table E1. Seasonal Adjustment Factors for Log Transformed PbB Data from 1990 
through February, 1996. 


Year 1990 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date aed aa 
Jan 1, 90 3.043 | g060 060 | 046 046 


Jan 16 3.060 3.080 0.077 


3.081 3.119 0.090 


Feb 16 3.121 3.164 0.097 


3.272 3.204 0.089 


Mar 16 3.431 3.251 0.114 


3.127 3217 0.146 


Apr 16 3.174 3.208 0.098 


3.179 3.212 0.088 


May 16 3.135 3.251 -0.057 


3.494 3.316 0.180 


Jun 16 3.249 3.339 -0.259 


3.536 3.351 0.300 


Jul 16 3.217 3.312 -0.263 


3.282 3273 0.207 


Aug 16 3.242 3.253 -0.210 


3.174 3.237 0.195 


Sep 16 3.265 3.249 -0.210 


3.289 3.247 0.200 


Oct 16 3.307 3.223 -0.159 


3.168 3.173 0.091 


Nov 16 3.003 3.124 -0.025 


3.088 3.123 -0.006 


Dec 16 3.187 3.143 -0.008 





77 


Table E1 (continued) - Year 1991 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date a =a 
Jan 1, 91 3.161 | gi59 | oot 


Jan 16 3.203 3.171 0.052 


3.159 3.170 0.106 


Feb 16 3.146 3.178 0.151 


3.164 3.202 0.181 


Mar 16 3.216 3.241 0.195 


3.304 3.292 0.133 


Apr 16 3.296 3.301 0.050 


3.250 3.308 0.032 


May 16 3.461 3.316 -0.113 


3.156 3.285 0.187 


Jun 16 3.324 3.276 -0.222 


3.288 3.272 0.252 


Jul 16 3.258 3.249 -0.223 


3.421 3.247 0.215 


Aug 16 3.375 3.256 -0.218 


3.241 3.285 0.191 


Sep 16 3.289 3.187 -0.137 


3.048 3.107 0.304 


Oct 16 3.002 3.030 -0.191 


2.97 2.973 0.098 


Nov 16 2.881 2.929 -0.017 


2.914 2.907 0.041 


Dec 16 2.858 2.897 0.087 
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Table E1 (continued) - Year 1992 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date ae aa 
Jan 1, 92 2.953 2. | a g97 fa aoa 2 101 


Jan 16 2.835 2.894 0.096 


2.914 2.899 0.082 


Feb 16 2.905 2.903 0.070 


2.951 2.896 0.068 


Mar 16 2.868 2.879 0.077 


2.809 2.876 0.072 


Apr 16 2.847 2.904 0.036 


2.924 2.964 0.032 


May 16 3.158 3.032 -0.107 


3.065 3.066 0.149 


Jun 16 3.098 3.079 -0.170 


3.029 3.088 0.186 


Jul 16 3.157 3.103 -0.209 


3.072 3.114 0.296 


Aug 16 3.146 3.120 -0.239 


3.169 3.114 0.241 


Sep 16 3.135 3.067 -0.201 


2.934 2.996 0.137 


Oct 16 2.968 2.934 -0.081 


2.817 2.875 0.028 


Nov 16 2.790 2.846 -0.006 


2.845 2.846 0.013 


Dec 16 2.865 2.841 -0.014 
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Table E1 (continued) - Year 1993 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date a a 
Jan 1, 93 2.923 2. | ago | 000 000 


Jan 16 2.787 2.757 0.058 


2.595 2.679 0.129 


Feb 16 2.601 2.635 0.168 


2.560 2.625 0.172 


Mar 16 2.637 2.650 0.141 


2.702 2701 0.085 


Apr 16 2.806 2.745 0.035 


2.728 2779 0.004 


May 16 2.873 2.820 -0.050 


2.784 2.852 -0.087 


Jun 16 2.928 2.897 -0.137 


2.969 2.935 0.180 


Jul 16 2.940 2.959 -0.209 


3.005 2.987 0.242 


Aug 16 2.958 2.993 -0.252 


3.091 2.996 0.260 


Sep 16 3.068 2.956 -0.225 


2.737 2.878 0.151 


Oct 16 2.839 2.836 -0.113 


2.804 2.805 0.087 


Nov 16 2.769 2.765 -0.050 


2751 2.729 0.019 


Dec 16 2.806 2.672 0.036 





80 


Table E1 (continued) - Year 1994 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date a ai 
Jan 1, 94 2.431 | a se7 | ooatz 


Jan 16 2.472 2.545 0.154 


2.578 2.547 0.149 


Feb 16 2.609 2.545 0.149 


2.455 2.544 0.146 


Mar 16 2.576 2.569 0.117 


2.577 2.604 0.080 


Apr 16 2.618 2.652 0.029 


2.766 2.732 0.054 


May 16 2.779 2.796 -0.121 


2.789 2.871 0.199 


Jun 16 3.143 2.942 -0.272 


2.926 2.952 0.285 


Jul 16 2.993 2.946 -0.281 


2.814 2.937 0.275 


Aug 16 2.958 2.946 -0.285 


3.010 2.962 0.303 


Sep 16 3.045 2.939 -0.282 


2.868 2.882 0.207 


Oct 16 2.806 2.808 -0.155 


2.652 2.748 0.097 


Nov 16 2.751 2.712 -0.062 


2.652 2.685 0.037 


Dec 16 2.741 2.658 -0.011 
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Table E1 (continued) - Year 1995 


Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date eae <a 
Jan 1,95 2.563 2. |p gig | 027 027 


Jan 16 2.586 2.576 0.069 


2.525 2.548 0.096 


Feb 16 2.591 2.522 0.122 


2.346 2.499 0.144 


Mar 16 2.549 2.514 0.128 


2.571 2.536 0.105 


Apr 16 2.560 2.549 0.091 


2.485 2.583 0.058 


May 16 2.606 2.647 -0.007 


2742 2.732 0.092 


Jun 16 2.901 2.808 -0.168 


2.899 2.842 0.202 


Jul 16 2.837 2.839 -0.199 


2.812 2.827 0.187 


Aug 16 2.782 2.815 -0.174 


2.800 2.816 0.175 


Sep 16 2.900 2.820 -0.178 


2792 2.812 0.170 


Oct 16 2.757 2.807 -0.164 


2.747 2.820 0.176 


Nov 16 2.989 2.843 -0.198 


2.867 2.824 0.178 


Dec 16 2.723 2.774 -0.127 
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Table E1 (continued) - Year 1996 


Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
eral Date an a 


Jan I = anttoes | 96 2.720 | azo | 721 | 9.072 072 


Jan 16 2.738 2.651 -0.001 


2.596 2.556 0.096 
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Table E2. Seasonal Adjustment Factors for 1990 to February, 1996 Based on Analysis of 
Untransformed Data. 


Year 1990 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date aaa al 
Jan 1, 90 | 2sc85 | 0.062 _| 


Jan 16 24.518 0.341 


sae | 0.256 


Feb 16 26.395 0.091 


27.350 | -0.060 


Mar 16 27.800 0.316 


zecoo | 1679 


Apr 16 gees | 1.173 


2ao25 | _o.7e7 


May 16 poss | 2.631 


oso | 4.529 


Jun 16 31.760 5.926 


31965 | 6.487 


Jul 16 30.965 5.712 


a760_ | 4.793 


Aug 16 pears | 367 


a7s10__ | 3.398 


Sep 16 ger | 4.10 


pases | 4.376 


Oct 16 27.355 3.285 


25900 | 1.967 


Nov 16 25.150 1.314 


2500 | 1.331 


Dec 16 25.145 1.541 
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Table E2 (continued) - Year 1991 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date eal al 
Jan 1,91 | aso10 | 1.108 _| 


Jan 16 24.570 0.044 


ago | 1.197 


Feb 16 24.010 2.031 


paaso | 2.275 


Mar 16 25.070 2.400 


e200 | 1273 


Apr 16 27.105 | 0.343 


paso | 2.267 


May 18 asso | 4.206 


esse | aig 


Jun 16 29.240 5.310 


pasio | 6.728 


Jul 16 29.130 4.937 


e960 | 4.354 


Aug 16 asss | 4.336 


ous | ari 


Sep 16 7845 | 1.998 


soo | 7am 


Oct 16 23.100 A711 


21.190 | 2.520 


Nov 16 19.395 0.443 


18.405 | 0.830 


Dec 16 18.045 1474 
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Table E2 (continued) - Year 1992 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date ae al 
Jan 1,92 | issa0 | 1270 _| 


Jan 16 18.620 0.888 


19130 | 0277 


Feb 16 19.250 0.056 


19.025 | 0.182 


Mar 16 18.935 0.173 


18.650 | 0.361 


Apr 16 1800 | 0.114 


20100 | 1.282 


May 18 ai4so | 2.727 


pono | 3572 


Jun 16 28.000 4.465 


paso | 6.007 


Jul 16 28.750 5.399 


4100 | 6.840 


Aug 16 24105 | 6.995 


paso | 8.729 


Sep 16 pao20 | 6.027 


21.4 21.930 -4.024 


Oct 16 20.820 3.001 


19560 | 1.926 


Nov 16 18.755 “1.108 


18200 | -0.635 


Deo 16 17.650 0.168 
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Table E2 (continued) - Year 1993 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date ae al 
Jan 1,93 | 17400 | 0.000 _| 


Jan 16 16.800 0.519 


15950 | 1.288 


Feb 16 15.450 1.709 


15.00 | _1.780 


Mar 16 15.700 1.302 


16370 | 0586 


Apr 16 e840 | 0010 


7630 | -0.856 


May 18 18470 | 1.770 


18730 | 2.103 


Jun 16 19.140 2.586 


9620 | 3.138 


Jul 16 20.150 3.738 


21.000 | 4.658 


Aug 16 i610 | 5.338 


pie0 | 5.486 


Sep 16 oss | 4.993 


1ac40 | 2.870 


Oct 16 18.125 2.120 


17590 | -1.660 


Nov 16 16.945 1.069 


16730 | -0.917 


Deo 16 16,280 0.529 
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Table E2 (continued) - Year 1994 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date el al 
Jan 1,94 | 14960 | 0730 _| 


Jan 16 14.245 1.384 


14240 | 1.990 


Feb 16 14.215 1.296 


14150 | 1.808 


Mar 16 14.450 0.946 


5410 | 0.230 


Apr 16 16170 | -0.985 


17790 | 2.569 


May 18 13945 | 3.768 


po47o | 6.348 


Jun 16 22.130 7.087 


pa2so | 7.238 


Jul 16 22.050 7.078 


aro | 6.767 


Aug 16 21555 | 6.681 


poo | 7.188 


Sep 16 a1es0 | 7.069 


posoo | 6.768 


Oct 16 18.650 3.960 


17950 | 2.704 


Nov 16 16.200 1.597 


15560 | -0.989 


Deo 16 15.350 0.830 
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Table E2 (continued) - Year 1995 
Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
Starting Date ae aaa 
Jan 1, 95 | iass0 | 0071 | 


Jan 16 13.900 0.539 


13560 | 0.861 


Feb 16 12.880 1.483 


12660 | 1.668 


Mar 16 12.970 1.320 


13.030 | 1.205 


Apr 16 13070 | 1.150 


13460 | 0.727 


May 18 14580 | -0.426 


16250 | 2.128 


Jun 16 18.095 4.004 


18790 | 4.729 


Jul 16 18.480 -4.448 


18925 | 4.901 


Aug 16 7940 | 3.964 


1110 | -4.160 


Sep 16 18.075 | 4.161 


7470 | 3571 


Oct 16 17.210 3.835 


7145 | 3.299 


Nov 16 17.330 -3.500 


7195 | -3:987 


Deo 16 16.810 3.022 
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Table E2 (continued) - Year 1996 


Time Period 90th Percentile PoB | Smoothed 90th | Adjustment 
il Date ae al 


| ant96 | tets0 | 282 | 


Jan 16 14.940 1.191 


13480 | 0282 
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Table E3. Age Adjustment Factors for 1990-February, 1996 Milwaukee PbB Data 


Adjustment for Log (PbB) Multiplicative Adjustment for 
Age Category (Years Values Untransformed Data 


0.5 - 0.75 0.812 2.252 


0.75 - 1 0.279 1.322 


1.0 - 1.25 0.096 1.101 


1.25 - 1.5 0.039 1.040 


1.5 - 1.75 -0.028 0.972 


175-2 fd 1.000 


25-3 1.040 
ee ee 0 
er es Ye 
ee eee 
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